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Abstract 

Introduction. A new scheme of a flat statically determinate regular lattice is proposed. The lattice rods are hinged. The 
study aims at deriving a formula for the dependence on the number of panels of the first natural oscillation frequency of 
nodes endowed with masses, each of which has two degrees of freedom in the lattice plane. The rigidity of all rods is 
assumed to be the same, the supports (movable and fixed hinges) — nondeformable. Another objective of the study is to 
find the dependence of the stresses in the most compressed and stretched rods on the number of panels in an analytical 
form. 

Materials and Methods. The approximate Dunkerley’s method was used to determine the lower bound for the lattice 
natural frequency. The lattice rigidity was found in analytical form according to Maxwell-Mohr formula. The rod 
stresses and the reactions of the supports were determined from the equilibrium equations compiled for all lattice nodes. 
Generalization of the result to an arbitrary number of panels was performed by induction using Maple symbolic math 
operators for analytical solutions to a number of problems for lattices with different number of panels. 

Results. The lower analytical estimate of the first oscillation frequency was in good agreement with the numerical solu- 
tion for the minimum frequency of the oscillation spectrum of the structure. Formulas were found for the stresses in 
four most compressed and stretched rods and their linear asymptotics. All required transformations were made in the 
system of Maple symbolic math. 

Discussion and Conclusions. The obtained dependence of the first frequency of lattice oscillations on the number of 
panels, mass and dimensions of the structure has a compact form and can be used as a test problem for numerical solu- 


tions and optimization of the structure. 
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Introduction. Lattice structures are widely used in mechanical engineering as load-bearing or enclosing elements. 
The procedure of calculating the natural frequencies of lattice vibrations is traditionally based on numerical evaluations, 
which, as a rule, use the finite element technique with application-specific packages. Analytical methods are used much 
less frequently. The development of mathematical packages of symbolic mathematics (Maple, Wolfram Mathematica, 
Derive, Maxima, etc.) enables to search for such solutions for regular systems in which it is possible to take into ac- 
count the order of regularity (the number of periodic structures of the construct, specifically, the number of panels) in 
the solution. The study aims at finding analytical dependence of the lower boundary of the first part of the lattice struc- 
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ture on the number of panels. The dependence of the solution on the number of panels expands significantly the scope 
of the formula and makes it possible to optimize the design by any parameter. For the first time, R. Hutchinson and 
N. Fleck were engaged in the problem of the existence and calculation of statically definable regular truss systems [1, 
2]. In [3, 4], such constructions were studied in connection with optimization problems. There are also analytical solu- 
tions in the form of finite formulas for the programs of regular planar [5—7] and spatial trusses [8]. Reference book [9] 
provides formulas for deflections and displacements of supports for flat trusses and lattices with an arbitrary number of 
panels. In [10—13], lower estimates of the first natural frequency of flat trusses were obtained, found by induction for an 
arbitrary order of a regular structure. 

There is another line of analytical research of structures [14—16]. In these works, the solution is in the Maple system 
in the form of trigonometric series. 

In this paper, the induction method has been used, which consists in generalizing a number of separate calculations 
of lattices with a consistently increasing number of panels for the case of an arbitrary number of panels. All transfor- 
mations used Maple symbolic math operators. The object of the study was a new scheme of a regular statically defina- 
ble flat lattice of rectangular shape in the form of an asymmetric truss on two supports. The task was to deduce the ana- 
lytical dependence of the fundamental vibration frequency of the structure on the number of panels. The derived formu- 
la can be used in optimization problems and for estimating numerical solutions of such large-order structures, for which 
numerical calculations may contain errors associated with the round-off accumulation. 


Materials and Methods. There are 5(7+1) nodes, including three support nodes (one for fixing the left support, 
two — for the right bars modeling a fixed hinge), in a truss with length (2n—Da and height 4h. The number of rods, 


including three support rods, is N =10n+4. The truss is statically definable. When calculating the vibration frequen- 


cies of the structure, it is assumed that the mass of the truss is concentrated in the nodes. 


Fig. 2. Lattice numbering, n=5 (the author's figure) 


The force calculation is performed in the Maple system according to the program [17]. The nodes and rods of the 
truss are numbered (Fig. 2). The origin of coordinates is in the left support. Coordinates are set in cycles. 
The lattice structure is established by the order of connection of the rods in the nodes. To do this, special lists 


®, =[i,,i,] of numbers i,,i, of the rod ends a=1,..,N are introduced. The rods of the lower outer contour, specifical- 
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ly, have the following node numbers at the ends: ®, =[i,i+1], i=1,..,n. The numbers of the ends of the other rods of 


the lattice are set in the same way. The system of equilibrium equations of nodes in projections on the coordinate axis is 
compiled in matrix form GS=B, S — vector of all rod stresses, including three reactions of the supports, B — vector 
of external nodal loads. System matrix G consists of guiding cosines of stresses. In this case, the same stress is applied 
to different rod ends and are multidirectional: 

Gyo, 14 = “hs /I,, Gyo, = —l ‘I; 

Gyo, = bi Ns Goo, =H il, T= LN, 


where 1. ; =X, —%o,,> 4); =Yo,, —Yo,, — Projections of conditional vectors of rods on the coordinate axis, 


Laat, +E, — length of rod i=1,..,N. 


The lattice rod stresses can be obtained from solving a system of equations in symbolic or numerical form. 

Research Results. Consider the stress state of the lattice in the case of loading at all angles by vertical 
stresses P (Fig. 1). Figure 3 shows a distribution pattern in the rods of the structure. The thicknesses of the segments in 
the Figure are conditionally proportional to the modules of the corresponding stresses. Compressed rods are highlighted 
in blue, stretched rods — in red. The stress values are attributed to magnitude P of the load on the node, rounded to two 
significant digits. The most stretched rod was expected to be in the middle of the lower belt, the most compressed 


one — in the lower rod on the right side of the structure. 


Fig. 3. Stress distribution on the lattice rods, a = 4m, h = 3m,n =5 (the author's figure) 


The analytical dependences of the stresses in the most compressed and stretched rods on the number of panels are 


obtained by induction from the generalization of the sequences of individual solutions. Say, for stress V,, in the lower 


ln 
rod of the left side of the lattice, the sequence of values for lattices of order n = 1, 2, 3, ... has the form: 


V,,/P=-2, —-12, —28, —54, —84, 126, —170, —228, —286, —360,.... The recurrent equation for the common 


term of this sequence is given by operator rgf_findrecur of the Maple system: 
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ln 


=Vi4 2V, 1-2 2V, 3 V, 


ln 


—4 Vins: 


Solving the equation using rsolve operator: 


V. =—P(30n? —2(7—(-1)")n—(-1)" +1) /(8(2n-D) . al7 
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Other expressions for critical stresses are found in the same way: 
O,= —Pa(20n* —30n? +5—4(cos @+sin @)(1— 2n) —cos(2@)) / (32(2n—- 1h); 
U, = =O 
V,, = P(5n? —n-1)/(2n-1), 


where @=7n/ 2. It is of interest to note that stresses V,,, and V,,, do not depend on the size of a and h. 


The analytical form of solutions allows us to find their asymptotics using Maple operator limit: 
limV,,, =—3P,,,/8, limV, ,, =—3P.,,,, /8; 


Pane sum on sum 


limO,/n=—lmU, /n=—aP.,, / (16h), 


sum 
noo n>o 


where P. =5nP — total load on the lattice. For stresses V, 


sum In 


and V,, in the rods on the sides of the lattice, the as- 


ymptotes are horizontal, for stresses O, and U,, on the upper and lower belts, they are inclined. 


Natural frequency. Of the entire spectrum of natural vibration frequencies of a structure, the first, lowest frequency 
is the most important for assessing its dynamic behavior. Its value is included in most solutions to the problems of struc- 
tural dynamics. This value is required, among other things, to assess the seismic characteristics of the structure. The 
lower bound of the first frequency for regular structures in the form of a dependence on the number of panels can be 
obtained analytically. 

When determining the natural vibration frequencies of the structure, a simplified model of the inertial properties of 
the truss is adopted. It is assumed that the rods of the lattice have no mass, and the entire mass is distributed evenly 
across the nodes. Neglecting the movement of the supports, we get the total number of degrees of freedom equal to 


K =10n. Donkerley’s formula [13] for estimating the lower bound of the first frequency has the form: 


K 
5 =>,» (1) 
p=l 
where @, — partial frequencies of the structure. Partial frequencies of mass vibrations are determined from the equa- 
tion: 
my, +D,y, =0, p=1,....K. (3) 
Here, y, = y,(t) — coordinate of node p; ¥,, — acceleration; D, — rigidity, inverse compliance 6, =1/D, . 


Compliance can be calculated using the Maxwell-Mohr formula: 


8, =1/D, = Y (sy 1, / (EF), (4) 
a=l 


where S‘’’ — stress in the rod with number « from the action of the vertical unit force applied to node p, where the 
mass is located. The rigidity rate and partial frequency depend on the place where the mass is located. For harmonic 
vibrations y, =U, sin(@t+@) from (3), @, =, /D, /m follows. The substitution of this expression in (4) gives a for- 


mula for estimating the first frequency only by partial frequencies of mass vibrations: 


K 


K 
> do, m>'6, mA, f Ag). (5) 


p=l p=l 


The amounts for vibrations in the vertical A, , and horizontal A,,,, directions are separately allocated. Successive 


v 


calculation of the vibration frequencies of lattices of various orders shows that coefficient A, , in (5) has the form: 
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A, , = (14h)/ (EF); 
A,, = (44a? +753 + 250h*) / (QEFh’); 

A, =(24a° + 25c? +36h°) / (EFh’); 

A,, = 2(154a° +175c* +176h°)/ (7 EFh’); 
A,, = 2(1254a? +1125c? + 791h°) / (27EFh’). 


ee) 


We 


In general, for vertical partial frequencies: 
Any =(C\a? +C,c° +C,,h°)/ (EFh’ ). (6) 


For partial frequencies horizontally: 


es) 


A, =(Cya’ + Cyc’ +C,,h*) /(EFa"). (7) 


Using Maple operator rgf_findrecur, we obtain a homogeneous recurrent equation of the seventh order to determine 


the coefficient for a’: 


C. 


In 


=3C,,4 Qs IC, 3 Cin Cee 3C,,,-6 Cunt . 
The solution to the equation is given by rsolve operator: 


C, =(n* = 2n’ +61n? —6(3(-1)" +10)n+9(-1)” —9)/18. (8) 


Other coefficients are found in the same way: 
C, = 25n(n-1)/6, 


9 
C, = (136n? + 2(9(-1)” -14)n-9(-1)" -15) /(6(2n-1)). ”) 


Similarly: 
G= (164n* — (36cos(2@) + 142)n? + (24cos(2¢) — 24 sine + 24cos p+ 20)n— 
—3cos(2o) + 12sin@—12cos@+ 27) /(2n—1)/12; 
C; = 10n* —5n? +15n+8)/(2(2n-1); (10) 
C, = (36n* — (28 +8cos(2@))n* + (128 + 22.cos(2¢@) + 24cos p— 24sin p)n— 
—12cos@—9cos(2@) +109+12sin@) /(2n—-1)? /2, 


where p= 2n/2. 
As a result, the expression for the lower estimate of the first frequency follows from (5—7): 
wp =m((Ca’ + Cyc + Ch) / 1? +(Cya’ + Cs! + Csh*)/ a” )/ (EF). (11) 
Estimation of the solution error (11) is possible from comparison to the minimum frequency of the entire spectrum 
of lattice natural frequencies obtained numerically. The spectrum of a system with many degrees of freedom is found 


from solving the eigenvalue problem of the matrix. Differential equations of mass dynamics of a structure with a num- 
ber of degrees of freedom K are written in matrix form: 


m,Y¥+D,Y =0, (12) 
where D, — truss rigidity matrix, Y— mass displacement vector, I, — unity matrix. Let B, be the matrix, inverse of 
D,.. Multiplying (12) from the left by B,, gives the equation: 

mB, ¥ +1,¥ =0. (13) 


Connection Y =—@°Y is valid for harmonic vibrations with frequency w. Hence, from (13), the problem of eigen- 
values follows: B, Y =AY, where 4 =1/(w°m) — the eigenvalues of matrix B, . The elements of the compliance ma- 
trix are found from the Maxwell-Mohr formula. To calculate the eigenvalues of a matrix in the Maple system, Eigen- 
values operator from LinearAlgebra package is used. In general, the solution to this problem can be obtained only in 
numerical form. 

Example. The steel lattice has panel length a=2m, height h=1m, masses in nodes m=300kg, rod rigidity 
EF =0.8-10°kN . Figure 4 shows the dependences of the first frequency on the number of panels obtained numerically 
and analytically. 
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1 2 3 4 5 6 7 8 9 10 1 12 #13 14 #+15 #16 #17 #18 =n 


Fig. 4. Comparison of the first natural vibration frequency of lattice ©, and its lower analytical estimate ©, 
depending on the number of panels (the author's figure) 


Relative error €=(®,—@,)/@,, depending on the number of panels (Fig. 5), shows that with an increase in the 


number of panels, the error, starting from m =1, increases, and then monotonously and quickly decreases. This is par- 
ticularly important when using the obtained formula through calculating lattices with a large number of panels, for 
which the numerical count starts to gain an error associated with the accumulation of rounding errors, and the cost of 


computer resources is growing rapidly. 
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Fig. 5. The error of the analytical solution for the natural frequency 
320 depending on the number of panels and height h (the author's figure) 
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Discussion and Conclusions. A scheme of a statically definable plane lattice was proposed. The task was to obtain 
an analytical expression of the lower bound of the first natural vibration frequency of a lattice truss for an arbitrary 
number of panels in the structure. The solution was obtained by induction in the Maple system. The presence of extrema 
on the constructed curves made it possible to optimize the number of lattice panels by choosing the highest accuracy of 
the estimate and selecting the required vibration frequency. The proposed study took into account horizontal mass vi- 
brations. Considering horizontal vibrations somewhat complicated the final calculation formula, making it more cum- 
bersome. In addition to the Donkerley’s method used, there was a more accurate Rayleigh energy method for evaluating 
the first frequency, which gave an estimate of the first frequency from above. However, even this solution in the case 


under consideration had an unnecessarily cumbersome appearance and was not given here. 
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